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1 Introduction 

Weak-scale supersymmetry (SUSY) [1, 2] predicts new particles with masses in the TeV 
range, which could be detected at the Large Hadron Collider (LHC). In particular the 
coloured squarks (g) and gluinos (g) would be produced abundantly in hadronic collisions. 
Searches for squarks and gluinos at the LHC have placed lower limits on squark and gluino 
masses around 1 TeV [3-6]. Once the LHC reaches its design energy near v5 = 14 TeV, 
SUSY particles with masses up to 3 TeV can be probed [7, 8]. 

In the context of the Minimal Supersymmetric Standard Model (MSSM) [9, 10] with im- 
parity conservation, squarks and gluinos are produced in pairs in collisions of two hadrons 
hi and /i2: 

hih2 -^ qq,qq,qg,gg + X . (1.1) 

In Eq. (1.1) and throughout the rest of this paper we suppress the chiralities of the squarks 
q = {qLiQfi) and do not explicitly state the charge-conjugated processes. We include 
squarks q of any flavour except for top squarks. The production of top squarks [11] has 
to be considered separately since the strong Yukawa coupling between top quarks, top 
squarks and Higgs fields gives rise to potentially large mixing effects and mass splitting [12]. 
Accurate theoretical predictions for inclusive squark and gluino cross sections are needed 
to set exclusion limits and to help determining SUSY particle masses and properties in 
case of discovery [13, 14]. The precision of the predictions can be improved considerably 
by including higher-order SUSY-QCD corrections, which have been known for quite some 
time at next-to-leading order (NLO) in SUSY-QCD [15-17]. A significant part of these 



corrections can be attributed to the threshold region, where the partonic centre-of-mass 
(CM) energy is close to the kinematic production threshold. The NLO corrections are then 
dominated by soft-gluon emission off the coloured particles in the initial and final state 
and by the Coulomb corrections due to the exchange of gluons between the slowly moving 
massive particles in the final state. The dominant contributions due to soft-gluon emission 
have the general form 

a^log'^/32 , m<2n with /32 = l-p=l ^, (1.2) 

where Og is the strong coupling, s is the partonic CM energy squared and rriav is the 
average mass of the final-state particles. These soft-gluon corrections can be taken into 
account to all orders in perturbation theory by means of threshold resummation techniques 
[18-23]. The all-order summation of such logarithmic terms is a consequence of the near- 
threshold factorization of the cross section. We perform the resummation after taking a 
Mellin transform (indicated by a tilde) of the hadronic cross section ah^h2^ki'- 

^hih2'^ki{N,{m^}) = dp/O^"^ o-/,,h2->fc/(/5,{?7T-^}) , (1-3) 

Jo 

with k and / the final-state particles and {m^} the squared masses involved in the process. 
The logarithmically enhanced terms are then of the form a^ log™ N, m < 2n, with the 
threshold limit /3 — t- corresponding to A^ — ?• 00. 

Threshold resummation has been performed for all MSSM squark and gluino pro- 
duction processes at next-to-leading-logarithmic (NLL) accuracy [24-28]. For squark- 
antisquark and gluino-pair production, in addition to soft-gluon resummation, the Coulomb 
corrections have been resummed both by using a Sommerfeld factor [25] and by employing 
the framework of soft-collinear effective field theories [29-31]. Furthermore, the domi- 
nant next-to-next-to-leading order (NNLO) corrections, including those coming from the 
resummed cross section at next-to-next-to-leading-logarithmic (NNLL) level, have been 
calculated for squark- antisquark and gluino pair-production [32, 33]. Very recently, ap- 
proximate NNLO predictions for stop-antistop production have been obtained using the 
soft-collinear effective field theory formalism [34]. For squark-antisquark production, soft- 
gluon emissions have been resummed to NNLL accuracy [35], resulting in a notable sta- 
bilization of the theoretical predictions. Recently, the NNLL resummation has been also 
performed for gluino-pair production [36]^. One would expect similar results for the other 
SUSY-QCD production channels. In the following, we discuss the calculation of a partic- 
ular class of ingredients necessary for performing NNLL resummation, the so-called hard 
matching coefficients. 

The squark and gluino production processes listed in Eq. (1.1) are scattering processes 
with a non-trivial colour structure. Since the soft radiation is coherently sensitive to the 
colour structure of the underlying hard scattering, colour correlations have to be taken 
into account when considering resummation. In an appropriately chosen colour basis, the 



^However, the result in [36] relies on results of [37] which we comment on in section 4. 



NNLL resummed partonic cross section f,-^^;.; takes the form of a sum over irreducible 
representations I [20, 29]: 

X A,(iV + l,Q2,/i2)A,-(iV + l,Q2,/i2)Ag)_^;-;,(Q/(iV^),/i2) . 

Here cf^^_^j^i j is the colour-decomposed leading order (LO) partonic cross section in Mellin- 
moment space, fi is the common factorization and renormalization scale, and we introduced 
the hard scale Q^ = 4?tt,^^. The last line of Eq. (1.4) captures all dependence on the large 
logarithm. The functions Aj and Aj sum the effects of the (soft-)collinear radiation from 
the incoming partons, while the function A- •\,^j j describes the wide-angle soft radiation. 
These functions are known at NNLL, see e.g. [35] and references therein. The matching 
coefficients in the second line contain the Mellin moments of the higher-order contributions 
without the log(A^) terms. This non-logarithmic part of the higher-order cross section near 
threshold factorizes into a part that contains the leading Coulomb correction C- °^^'[y and 

a part that contains the NLO hard matching coefficients Q_^fcj/ [20, 29]. These hard 
matching coefficients follow from the threshold limit of the full NLO calculation. They are 
a key element in the NNLL calculation and will be the subject of this paper. 

We will present the colour-decomposed NLO hard matching coefficients for the squark 
and gluino production processes at the LHC and discuss the impact of the corrections. 
We will start by constructing an appropriate colour basis in SU(A''c) in section 2 and 
discuss why some colour-decomposed cross sections vanish at threshold in section 3. In 
section 4 we discuss the calculation of the matching coefficients. The numerical results are 
presented in section 5 where we show predictions for the LHC with centre-of-mass energy 
of v5 = 8 TeV. We conclude in section 6. 

2 Construction of the colour bases in SU(Nc) 

When performing resummation for coloured particles, one has to take into account the 
colour correlations introduced by gluon radiation. In particular, wide-angle soft radiation 
is coherently sensitive to the colour structure of the hard process from which it is emitted 
[20-23, 38, 39]. At threshold, the resulting colour matrices become diagonal to all orders 
by performing the calculation in an s-channel colour basis [24, 25, 40]. This basis traces 
the colour flow through the s-channel and is obtained by performing an s-channel colour 
decomposition of the reducible two-particle product representations into irreducible ones. 
Methods to obtain such a basis have been presented in [26] and [40] for processes that 
contain particles of the same representations in the initial and final state and recently for 
any number of partons in SU(A^c) in Ref. [41]. Here we will present a method that is also 
valid for processes that contain particles of different representations in the initial and final 
state. All results are shown for a general SU(A'^c) theory with Nc the number of colours. 



We denote the colour charge operator of a representation by Tj. It is given by gluon 
emission off the corresponding particle j, so 

{'^q,g)ab = ^afe) i'^q,q)ab = ~'^ba = C^abY ^^'^ {Tg,g)ab = ^ab = —^fabc , (2.1) 

with T^j and F^j^ the generators of the fundamental and the adjoint representation respec- 
tively. In addition to the completely antisymmetric structures F^^, we will also need the 
traceless symmetric octet structures D^^ = dabc and the singlet colour structures 6ab- The 
dimension of the colour labels is determined by the particle they refer to. In SU(A'^c)i 
the quarks and squarks are in the Nc-dimensional fundamental representation, while the 
gluons and gluinos are in the adjoint representation with dimension N^ — 1. The colour 
decompositions for the squark and gluino production processes are given by: 

qq^q^: 1 © (N^ - 1) , (2.2) 

gg^m- lffi(N2-l)A©(N2-l)s, (2.3) 

qq^M- le(N2-l)A©(N2-l)s, (2.4) 



gg^gg- 1 (N^ - 1)^ (N^ - l)s e (N^ - 1){nI - 4)/4 (N^ - 1)(N2 - 4)/4 
N2(Nc + 3)(Ne - l)/4 e N2(Ne - 3)(Ne + l)/4 , (2.5) 

qq^qq- Ne(Ne - l)/2 © Ne(Ne + l)/2 , (2.6) 



qg^qg: Ne Nc(Ne + l)(Nc - 2)/2 © Ne(Ne - l)(Ne + 2)/2 . (2.7) 

In SU(3), the Nc{Nc — l)/2-dimensional representation for the qq — )• qq process coincides 
with the antifundamental representation 3, while the last representation of the gg — )• gg 
process vanishes. 

An s-channel colour basis can be constructed from any complete basis by requiring 
orthogonality, normalization and that it is an eigenvector under the quadratic Casimir 
operator. The orthogonality and normalization follow from an inner product of two colour 
tensors c/ and cj that describe the colour content of the 2 — )• 2 process under consideration: 



ci ■ cj = (^'J C^j) = 0/(01,02,03,04) c}(ai, 02, 03,04) = dim(c/)5/j, (2.^ 




where oi and 02 are the colour labels of the initial-state particles and 03 and 04 the colour 
labels of the final-state particles. The last equality in Eq. (2.8) fixes the normalization as 
well as the orthogonality. In addition, base tensors have to be eigenvectors of the quadratic 
Casimir operator (Tj + Tj)^, with Tj and Tj the colour charge operators of the two initial- 
state particles: 

((^i)aid/a2<i2 + (.'^j)a2d2^aidij (^(^i)dib/rf2b2 + (^j)d2fe2'^'iifei j '^a3fe3'^a4b4C/(bl, ^2, ^3, ^4) 

= C2{Ri)ci{ai, 02, 03, 04) , 



or in a shorter notation: 

{T, + Tjf®ci = C2{Ri)ci. (2.9) 

Here C2{Ri) is the quadratic Casimir invariant of the representation Ri that corresponds 
to the base tensor c/. Eq. (2.9) provides both the value of the quadratic Casimir invariant 
and an additional constraint for the base tensor. 

The combined requirements of orthogonality, normalization and invariance under the 
quadratic Casimir operator fix the base tensors up to a phase. The resulting base ten- 
sors and their corresponding dimension and quadratic Casimir invariant for SUSY-QCD 
processes are listed in Appendix A. 

3 Suppressed cross sections 

Near threshold, the s-wave contribution to the final state dominates. Higher values of 
the final-state orbital angular momentum quantum number Lgn are suppressed by higher 
powers of /?. Thus a cross section can be regarded as being "suppressed" near threshold if 
the Lfin = mode is not accessible due to symmetry properties. In this section, we will first 
list these general symmetry considerations and then apply them to the different processes. 
If the initial or final state consists of two identical particles, the symmetry properties 
are determined by the eigenvalue P of the permutation operator that interchanges the two 
particles: 

P = S^ (_l\L+S-S^-S2 = / +^ ^^"^ identical bosons, 

\ — 1 for identical fermions. 

Here L is the total orbital angular momentum of the particle pair in the CM frame, S the 
total spin, 5i and 5*2 the spins of the individual particles and Sc the colour symmetry 
factor, which is +1 for a symmetric colour state and —1 for an antisymmetric colour 
state. Because we have Lgn = for a non-suppressed cross section, the conserved angular 
momentum quantum number equals the final-state spin, i.e. J = S^^. In addition, we 
will use the conserved total angular momentum quantum number Mj G { — J, . . . , J} for 
quantization along the propagation axis of the incoming massless partons (z-axis). In that 
case Mj is given by the difference of the helicities of the initial-state partons. 

QiQi -^ QiQi- only the symmetric colour structure contributes. We will argue that 
if the produced squarks in the qq — t- qq process have the same flavour, the contribution 
from the antisymmetric colour structure is suppressed near threshold. In that case, we 
have a system with identical fermions in the initial state with P = —1 and identical bosons 
in the final state with P = +1. Also, since squarks are scalars, the total spin of the 
final state is Sfin = 0. We need Lfin = for a non-suppressed threshold cross section, so 
the conserved angular momentum quantum number is J = 0, which automatically means 
Lin = Sin- Inserting this into Eq. (3.1) for the initial-state quarks, we see that the exponent 
is always odd. Eq. (3.1) shows that we can only have an antisymmetric state if 5c = +1) 
i.e. the colour structure is symmetric. Indeed, as we set out to argue, for equal flavours 
the antisymmetric colour structure is suppressed near threshold. 



This latter statement can be derived even more directly from the final state, which 
should be symmetric since squarks are spin-0 particles. Given that 5fin = 5i = 52 = and 
that we want Lg^ = 0, one readily obtains the requirement that Sc = +1 from Eq. (3.1). 

gg -^ qq: only the symmetric colour structures contribute. In a similar way as 
in the first argument for qiQi — t- qiQi, we can show that for the gg — t- qq process only the 
symmetric colour structures contribute near threshold. In this case the gluons need to be 
in a symmetric state, since they are spin-1 particles. As with the previous process, we have 
J = based on the final state. Using Eq. (3.1) we can see that the exponent is always 
even, so also in this case only the symmetric colour states can yield a non-suppressed 
contribution near threshold. 

gg — )■ gg: only the symmetric colour structures contribute. First, note that for 
the production of a pair of gluinos the final state should be antisymmetric, since gluinos are 
spin-i particles. Given the fact that we need Lgn = (see above) and that J = Sa^ = or 
1, Eq. (3.1) leads to the requirement that {—l)'^~^Sc = —1. So, if the gluinos are produced 
in a J = (J = 1) state the colour structure should have Sc = -|-1(— 1), i.e. be symmetric 
(antisymmetric). Which of these situations is realized will depend on the initial state. 
For the gg — )• gg process, the gluons have helicities ±1 and the total angular momentum 
quantum number for quantization along the z-axis is accordingly given by Mj = ±2, 0. 
The case Mj = corresponds to gluons with equal helicities and the case Mj = ±2 implies 
J > 2, which is ruled out by the final-state requirement that J = or 1. This J-ambiguity 
we must now settle. 

At this point we would have liked to invoke the Landau- Yang theorem [42] to rule out 
the J = 1 option. After all, in the threshold limit the gluino pair is produced at rest 
and can therefore effectively be regarded as a spin-J particle that (in the time-reversed 
sense) is decaying into two massless spin-1 gluons. However, if we carefully follow the 
argumentation of the Landau- Yang theorem, which is formulated for photons rather than 
gluons, no additional constraint is obtained. The colour quantum numbers of the gluons 
provide in fact a loophole for evading the implications of the Landau- Yang theorem. 

Nevertheless in our results we observe that only the symmetric colour states yield a non- 
suppressed threshold contribution to the LO and NLO cross sections, implying that only 
the J = state survives. The reason for this turns out to reside in the LO matrix element, 
which also features prominently in the calculation of the NLO virtual corrections. This 
LO matrix element can be cast into a multiplicative form consisting of two factors, one 
containing the colour structure and one containing all Lorentz and spin structure. In view 
of the gluino-pair symmetry arguments presented above, this means that only symmetric 
or antisymmetric colour structures can survive in the threshold limit and not both types 
of structures simultaneously. 

qq -^ gg: only the antisymmetric colour structure contributes. For the qq — t- gg 

process we can make use of the fact that the helicities of the massless initial-state quark 
and antiquark are opposite, as a result of chirality conservation in SUSY-QCD. The case of 
equal helicities is suppressed by the negligible mass of the initial-state quarks. Therefore, 



Mj = ±1 for angular momentum quantization along the z-axis, which excludes the J = 
gluino-pair state. Based on the gluino-pair symmetry arguments presented earlier, this 
requires the gluinos to be in an 5c = —1 state, i.e. an antisymmetric colour state, as is 
indeed observed. 

These symmetry arguments remain also valid at higher orders, so the hard matching 
coefficients of suppressed cross sections only contribute to 1/A^-suppressed terms. ^ Since 
these higher-order 1/A^-suppressed effects are beyond the level of accuracy of resummation 
we want to perform, we put the value of the hard matching coefficients for the suppressed 
cross sections to zero. Note however, that we choose to keep the 1/A^-suppressed terms 
originating from the LO cross sections in Eq. (1.4). 

4 Calculation of the hard matching coefficients 

In this section we will discuss the calculation of the hard matching coefficients C^^' at one 
loop. The hard matching coefficients C'^^ are determined by the terms in the /3-expansion 
of the NLO cross section that are proportional to /3, /31og(/3) and /31og^(/3). These terms 
receive contributions from both the real and the virtual corrections. After taking a Mellin 
transform, terms that contain higher powers of /3 are suppressed by powers of l/yiV and 
do not have to be considered. 

To obtain the integrated real corrections at threshold, we observe that they are formally 
phase-space suppressed near threshold unless the integrand compensates this suppression. 
Therefore we can construct the real corrections at threshold from the singular behaviour of 
the matrix element squared, which can be obtained using dipole subtraction [43, 44]. This 
has been discussed in detail in Ref. [35] in the context of squark-antisquark production. 
The same procedure can be used for the other SUSY-QCD processes and the result in 
n = 4 — 2e dimensions is given by: 

af '*"■■ = l67ra,Sna^''^'''4c2{Ri) (^ - logCS/?^) + 3) (4.1) 



n={i,j} 



^ ^ ^og(2) - -^ ) + log2(8/32) - 41og(-«2^ 



2^2 e I ^^ ^ 2r2 



+«-¥-'°«(lfc)0°^(«^^>-^+2&-'°^(^) 



"av 



Thus the real threshold cross section a^ ' ^ corresponding to a representation / is pro- 
portional to the colour-decomposed LO cross section at threshold aj ' ^ . The final-state 
contributions are weighted by the quadratic Casimir invariant C2{Ri) of the representation 
i?/, which is listed in Appendix A. The sum in the last two lines of Eq. (4.1) runs over 



^An exception might be the gg — >■ gg case, which might receive a contribution from the antisymmetric 
colour structures at NNLO level from squaring the NLO matrix element. However, the hard matching 
coefficients at NLO are not affected by this. 



the initial-state partons. The colour operators T^ depend on the representation of the 
corresponding particle, while the values of the flavour coefficients 7„ are determined by the 
partons in the initial state: 

where n/ = 5 is the number of light flavours. Finally, fxp is the factorization scale and the 
factor Sn is given by: 



^n 



167r2 



.2 



^ „-e(7B-log(47r)) / f^R 



m 



2 



"av , 

with ^E being Euler's constant and /i/j the renormalization scale. Eq. (4.1) is valid for all 
pair-production processes in SUSY-QCD that have a linear behaviour in /? at LO, including 
the unequal mass case of squark-gluino production. 

To obtain the virtual corrections, we start from the full analytic calculation as presented 
in Ref. [17]. As described in detail in Ref. [17], the QCD coupling a^ and the parton 
distribution functions at NLO are defined in the MS scheme with five active flavours, with 
a correction for the SUSY breaking in the MS scheme. The masses of squarks and gluinos 
are renormalized in the on-shell scheme, and the top quark and the SUSY particles are 
decoupled from the running of a^. 

We first need to colour-decompose the virtual corrections and then expand them in 
/3. These steps have been discussed for squark-antisquark production in Ref. [35] and the 
same procedure can be followed for the other processes although for the processes involving 
gluinos one does need to pay attention to spurious 1//3 singularities originating from Gram 
determinants occurring in the reduction of tensor integrals. After expanding to sufficient 
powers in /3, these singularities cancel in the calculation. The hard matching coefficients 
are now obtained by adding the non-logarithmic terms of the Mellin transform of Eq. (4.1) 
to the Mellin transform of the Oi^ji) virtual corrections to the cross section and dividing 
the result by the LO threshold cross section in Mellin-moment space. Note that the NLO 
hard matching coefficients will thus also contain non-logarithmic contributions from Mellin 
transforms of the 0(/3) logarithmic terms in Eq. (4.1). The complete expressions for the 
hard matching coefficients of the SUSY-QCD processes can be found in Appendix B. The 
virtual part of the corrections to the gg — )■ gg process agree with the results presented 
in [37], provided that one translates their DR result to our MS result and manually de- 
couples the heavy particles from the running of Og. On top of that, the top-quark mass 
has been neglected with respect to the squark and gluino masses in the results of Ref. [37] . 
However, the numerical impact of this approximation is small. Applying the same trans- 
lation procedure to the virtual part of the corrections to qq ^ gg., we find an agreement 
with the results presented in [37], apart from the expression for C4(r) in Eq. (32) of [37], 
where the ln( — ^ 1 term should read Inf 1 + r J .^ 



^We thank P. Marquard for confirming this typo in Eq. (32) of [37]. 



4.1 Numerical properties of the hard matching coefficients 

The behaviour of the hard matching coefficients for a varying squark-gluino mass ratio 
r = mg/rriq is shown in Fig. 1. The corresponding plots for squark-antisquark production 
can be found in Ref. [35]. 
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Figure 1. Mass dependence of the SU(3) colour-decomposed NLO hard matching coefficients for 
the qq — > gg process (a), the gg — >■ gg process (b), the qq — )■ qq process (c) and the qg — >■ qg process 
(d). The common renormalization and factorization scales have been set equal to the average mass 
of the produced particles vfiav ~ 1-2 TeV. The top quark mass is taken to be rrit = 172.9 GeV. 



We see from Fig. 1 that the hard matching coefficients can be quite large depending 
on the masses involved. The processes that involve gluinos have a singularity if the gluino 
mass equals the sum of the squark and the top mass. This singularity could be removed by 
taking the finite gluino width into account. These matching coefficients have been checked 
numerically using PROSPINO [17]. As explained in section 3, the hard matching coefficients 
of suppressed cross sections have been set to 0. 

As we can see from Fig. 1(a), the hard matching coefficient C ^^-- § for the qq — )• gg 
process is ill-defined at r = 1. This behaviour originates from the threshold behaviour 
of the LO cross section around r = 1. For the qq — )• gg process, the leading term in the 



large A'^ expansion is proportional to (m? 



m?)^ at LO and thus vanishes accidentally 



for VHn 



ms. 



The first terms in the large N expansion for the non-suppressed LO cross 



section in Mellin moment space are given by: 



1 + 2B+39B'^+34.B^-12B^-AB^+9B'^ 



^^aa sAN)^ "-" .'^'%/;' ( B^ - -^--^ B' (4.2) 



+ - 



4(iV+§)(iV+i 



with B = ^vTj which vanishes for r = 1. Thus the first two terms vanish for equal squark 
and gluino masses and the first non-zero term in the LO cross section is the 0{N~'^''^) 
term, which corresponds to a 0{I3^) threshold behaviour. There is no symmetry that 
causes this behaviour, so it is not surprising that it does not hold at higher orders. In 
fact, the NLO threshold cross section contains terms proportional to m? — tti? To obtain 
the hard matching coefficient from the NLO threshold cross section, we have to divide the 
NLO cross section by the 0{N~^''^) leading term in Eq. (4.2), resulting in a (?7i? — w,?)~^ 
divergence in the hard matching coefficient. 

Although the hard matching coefficient is multiplied by a cross section that is sup- 
pressed near r = 1, such behaviour leads to numerical instabilities. For a stable numerical 
implementation, we thus propose to define a modified matching coefficient, which uses 
higher order terms in 1/A^ to regularize the divergence. Instead of using only the leading 
term in the large A^ expansion, we use all three terms given in Eq. (4.2). After working 
out the Mellin transforms, the modified hard matching coefficient takes the form: 



^(l,mod) 



A+B-B"^ 1 + 2B+39B'^+3AB^-12B'^-AB^+9B^ 

AT+I ^ 4S2(Ar+|)(Ar+f) 






(4.3) 



This deeper expansion of the LO cross section ensures a well-behaved hard matching coef- 
ficient, which vanishes at r = 1. The hard matching coefficient now depends on A^, which 
is a complex variable. In order to obtain the full cross section, we integrate over an ap- 
propriate contour in A^-space, for which we use the "minimal prescription" of Ref. [45]. 
The behaviour of both the real and the imaginary part of the hard matching coefficient 
is shown in Fig. 2. For large values of |i?e(A^)|, the modified hard matching coefficient 
approaches the value of the A^-independent unmodified matching coefficient, cf. Eq. (4.3), 
and the curves get fiatter with increasing |i?e(A^)| in Fig. 2. As r approaches 1, the second 
term in Eq. (4.3) eventually takes over and Q^^L gA becomes (not shown in the plots), 
also for large (fixed) |i?e(A^)|. 

5 Numerical impact on LHC cross sections 

In this section we illustrate the numerical impact of the hard matching coefficients on 
the NLL-resummed cross sections. For this we first have to perform the inverse Mellin 
transform in order to recover the hadronic cross section crhih2~^kl- Iii order to retain the 
information contained in the complete NLO cross sections [17], the NLO and resummed 
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Figure 2. Real and imaginary part of the modified non- vanishing qq — > gg hard matching coefficient 
in tlic / = 8A colour channel for different values of r and TV = 2 -I- exp (i37r/4)C, C = O..00. The 
common rcnormalization and factorization scale has been set equal to the average mass of the 
produced particles niav = 1-2 TeV. The top quark mass is taken to be rrit = 172.9 GeV. 

results are combined through a matching procedure that avoids double counting of the 
NLO terms: 






(5.1) 



E 



CT 



2iri 



S:r"="''(iv.{™n.M^) - ^r:r"""''(A'.{™^).A.^) 



I (NLO) 



Here fi/h-^ and fj/h2 ^^^ ^^^ Mellin transforms of the parton distribution functions. We 
adopt the "minimal prescription" of Ref. [45] for the contour CT of the inverse Mellin 
transform in Eq. (5.1). 

We show the results for the LHC for a centre-of-mass energy of 8 TeV. In order to 
evaluate hadronic cross sections we use the 2008 NLO MSTW parton distribution func- 
tions [47] with the corresponding as(M^) = 0.120. We have used a top quark mass of 
mt = 172.9 GeV [48]. The numerical results have been obtained with two independent 
computer codes. In order to evaluate resummed cross sections, we use both the method 
of Ref. [49] with standard parametrization of pdfs in (c-space, as well as the method first 
introduced in [27], relying on Mellin-space pdf obtained with the program PEGASUS [50]. 
For comparison, we will also show results for the NLL matched cross section (jNLO-i-NLL 
without the effect of the hard matching coefficient, as well as LO and NLO results. 

The NLO cross sections are calculated using the publicly available PROSPINO code [51], 
based on the calculations presented in Ref. [17]. The QCD coupling Og and the parton 
distribution functions at NLO are defined in the MS scheme with five active flavours. 
The masses of squarks and gluinos are renormalized in the on-shell scheme, and both the 
top quark and the SUSY particles are decoupled from the running of ttg and the parton 
distribution functions. No top-squark final states are considered. We sum over squarks 
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with both chirahties {ql and gjj), which are taken as mass degenerate. The renormahzation 
and factorization scales /i are taken to be equal. 

We define the ET-factors for the NLL matched cross section with and without the hard 
matching coefficient in the following way: 



K^ 



a 



a(NLO) ■ 

These iT-factors are shown in Figure 3 for different mass ratios r = mg/rriq. We see that 
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Figure 3. The if-factor with respect to the NLO cross section of the NLO+NLL and NLO+NLL+C 
cross sections for the different SUSY-QCD processes for the LHC at 8 TeV. The common renor- 
mahzation and factorization scales have been set equal to the average mass of the pair-produced 
particles m = rriav 

the hard matching coefficients give a correction to the cross section that is comparable 
to the NLL correction. The corrections range from a few percent to up to 10% in the 
low-mass case, whereas they can reach several tens of percent in the high-mass case. Again 
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the largest effects are observed for those final states that give rise to the largest colour 
factors. In particular for processes involving gluinos, it is vital to use resummed results for 
obtaining an accurate prediction of the cross section, with total resummation corrections 
that can exceed 100% of the NLO cross section. The results are in general mildly mass- 
dependent, although a different mass ratio can lead to significant changes for gluino-pair 
production, and to a lesser extent for squark-pair production. The changes in the observed 
dependence on the mass ratio follow the r-dependence of the hard matching coefficients 
shown in Fig. 1. 

Although a full NNLL analysis is needed to draw definite conclusions on the size of the 
cross section and the scale dependence, these results indicate that the impact of the hard 
matching coefficients is quite significant. Thus we expect a fully NNLL-resummed cross 
section will improve the current theoretical predictions considerably. 

6 Conclusions 

In this paper, we have completed the calculation of the hard matching coefficients for 
the squark and gluino hadroproduction processes. The hard matching coefficients are an 
important ingredient for NNLL resummation. Numerically, we find that the inclusion of 
the hard matching coefficients increases the cross section at the central scale. 

Acknowledgments 

This work has been supported in part by the Helmholtz Alliance "Physics at the Teras- 
cale", the Foundation for Fundamental Research of Matter (FOM), program 104 "Theoret- 
ical Particle Physics in the Era of the LHC" , the DFG SFB/TR9 "Computational Particle 
Physics" , Polish National Science Centre grant, project number DEC-2011/01/B/ST2/03643, 
European Community's Marie-Curie Research Training Network under contract MRTN- 
CT-2006-035505 "Tools and Precision Calculations for Physics Discoveries at Colliders" 
and by the Research Executive Agency (REA) of the European Union under the Grant 
Agreement number PITN-GA-2010-264564 (LHCPhenoNet). 

A Base tensors for SUSY-QCD 

In this appendix, all the base tensors needed for the 2 — )• 2 SUSY-QCD processes are 
listed. They are obtained with the method described in section 2 and are given for a 
general SU(A''c) theory in terms of Kronecker deltas in colour space 6ab, the generators 
of the fundamental representation T^^, the structure constants fabc and their symmetric 
counterparts dabc- We give the dimension and quadratic Casimir invariant of all the base 
tensors. We use colour labels ai and 02 for the initial-state particles and labels 03 and 04 
for the final-state particles. Summation over repeated indices is implied. 

For the process q{ai)q{a2) — )■ 9(03)9(04) we have: 

Cqq^qll = J^Sa^a2^a3ai , dim(i2i) = 1 , C2(i?l) = , (A.l) 
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C,,-^,t,2 = 2 KaXsa, , dim(i22) =N^-l, C2{R2) = N, . (A.2) 

In SU(3), we have i?i = 1 and R2 = 8. 



For the process g{ai)g{a2) -^ Q{cLs)q{a4) there are three structures: 
1 

V^c(ivi-i)- 



Cgg->g|,l = /.,...., ^^aia2^a3a4 > dim(i?i) = 1 , C2(i?l) = , (A. 3) 



Cgg-^^12 = ^Ja,a,cT^,a, , dim(i?2) = iV ^ - 1 , C2(i?2) = N, , (A.4) 



V2]V; 



Cgg-,m = J-^ da,a,cT^,a, , dim(i23) = iV,^ " 1 , C2(i?3) = iV, (A.5) 



where Ri = 1, R2 = 8a and i?3 = 8s m SU(3). 

For the squark-pair production process q{ai)q{a2) — ^ ^(03)^(04) the base tensors are: 

Cqq-^qq,l ^ 7T ("0.10.4 '^012013 ~ ^ aia-j^ a2a4,) 1 C.qq^qq^2 ^ 77 (.0010400203 + 0aia30a2O4J j l-^-OJ 

and their dimension and quadratic Casimir invariants are given by: 

dim(i?i) = iiV,(Ar^ - 1) , C2(i^i)= ^^^^^^^^^"^\ (A.7) 

dim(i?2) = \n.{N, + 1) , C2(i22) = ^^' ~ ^^^^^ ^ ^^ . (A.8) 



'c 



In SU(3), Ri = 3, while R2 = 6. 

The colour structures of the q{ai)q{a2) —5- 5(03)5(04) process are similar to Eqs. (A.3- 

A.5): 

Cqq^gg,! = ^======(^ai02 5a304 , dim(i?i) = 1 , C2(i?l) = , (A.9) 

V^c{N^ - 1) 

C,,-^^3,2 = 7^/o304cT,^2ai , dim(i?2) = iV^ - 1 , C2{R2) = N, , (A.IO) 



Cqq^gg-S = -^^=da,a,cTl^a, , dim(iZ3) = iV^ - 1, C2{Rz) = N, . (A.ll) 

In SU(3), we have Ri = 1, R2 = 8a and i?3 = 8s. 

For the 5(01)5(02) — )■ 5(03)5(04) process, the base tensors are given by: 

^99^39,'^ ~ t\t2 _ 1 Oaia2 0a3a4 j (^-1^) 

^99^99,"^ ~ t\t2 _ A"'aia2cdca3a4 , (A. 13) 
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^gg^^gg^'i — at- Jaia,2cJcasa4, ) {■^•^'^) 



'c 



— L (X A —A A "l— ■^"iQ2c/ca3a4 i^/'j f if j N 

^gg-^ggA — a 10^^^300204 0Q,ja4 00203 j „,x 1 ^ VOoia3c/co204"r/oi03cOca204J ) 

(A.15) 

- — — ( X A —A A "l— J "■jo-'jcj ca^ai _^/j f .f ^ \ 

gg^fgg^^ — a \"a\a30a2aA 0Q,ja40o2O3J „ , x . (,"01030/00204 "r/aia3c"co204J ) 

(A.16) 

Nc + 2 ^ , A^c + 2 ,^ A ^A A ^ rA 17^ 

^gg^ggfi ~ oat (at , i\"ai«2 0o3O4 + ,^ (,00^0-300204 + 0aiO40o2O3J (,^--'-'/' 

iVc + 4 1, . 

"Uoia2cUa3a4C + "7 (,Wai03c'^0204C ~r 0'a2a3cWaia4cj ; 



Nc — 2 Nc — 2 

^gg^gg,'^ ~ o^v^T/V^^^TT " "' 4/V v0aia30a2O4 + 0010400203) (A.loj 



AtT\T r)\ ^0i02c'^a3a4C i (,"aj^o3cUo204C I Oio2a3cOaj^a4cj ; 



C2{Ri) = ^, 




(A.19) 


C2{R2)=N,, 




(A.20) 


C2{Ri) = N,, 




(A.21) 


C2{Ra) = 2Nc, 




(A.22) 


C2{R5) = 2N,, 




(A.23) 


C2iRe) = 2{N, + l), 


(A.24) 


C2{R7) = 2{N, - 


-!)• 


(A.25) 



^c-4 ^ _1 

while the corresponding dimensions and quadratic Casiniir invariants are: 
dim(i2i) = 1, 
dim(i22) = Nl-l, 
dim(i?3) = Nl-l, 
dim(i?4) = (iVe'-4)(A^,'-l)/4, 
d\Tn{R,) = {Nl-A){Nl-l)/A, 

dim(i?6) = iVc (A'c + 3)(iVe - l)/4 , 
diYn{R^) = Nl{N, - 3){N, + l)/4 , 

Note that since the dimension of Rj vanishes for Nc = 3, this representation does not 
contribute in (SUSY)-QCD. The other representations correspond to i?i = 1, -R2 = Sg and 
R3 = 8a- The i?4 and R^ are the 10 and 10 representations, while Rq = 27 in SU(3). 

Finally, the base tensors for squark-gluino production 9(01)5(02) — )■ q{a3)g(ai) are 
given by: 

2A^c 
^qg^qgA — j\^2 _ i^"^ ^"^ ^Jaaai ) (A. 26) 

Nc — 2 Nc — 2 

Cqg-^qg,2 = „ ^ 0020400403 ~ "'Ca4a2 -'^ a^ai + ]\[ _ -\ ^ Jaso-i ) (A. 27) 

Cqg-^qg,2 = „ /U "020400103 + "CO4O2 -'0301 ~ y^ _|_ 1 '^ )«3ai ' (A. 28) 

The dimensions and quadratic Casimir invariants for the corresponding representations 



are: 



dim{Ri) = Nc, C2(iii) = %-^, (A.29) 



2Nc 
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dim(i?2) = \Nc{N, + 1){N, - 2) , C2{R2) = (^^ ^l^^^"^^'^ , (A.30) 

dim(i?3) = \Nc{Nc - l)(iVc + 2) , C72(i?3) = ^^'^^2)^^'"^^ ' (A-31) 

which correspond to i?i = 3, i?2 = 6 and i^s = 15 in SU(3). 

B Hard matching coefRcients for SUSY-QCD 

Here we present the exact expressions for the hard matching coefficients C^^> for the SUSY- 
QCD production processes. We sum over squarks with both chirahties {qL and qn). No 
top-squark final states are considered and ah squarks are considered to be mass-degenerate 
with mass rriq. Top squarks are taken into account in the loops, where they are taken 
to be mass-degenerate with the other squarks. This is a valid approximation, since SUSY 
parameter dependence coming from loop corrections is generally small [17]. The calculation 
is outlined in section 4 and was done with FORM [52]. We first define the functions: 



t ■i\ /1 4mim2 , 2^ Mq ) - 1 



and the combinations: 



2 22 2 2,2 

m_ = rrig — m^, m_,_ = rrig + m^, 



where rrig is the gluino mass and rriq the squark mass. Furthermore, niav is the average 
mass of the produced particles. Denoting the number of light fiavours by n/ = 5, the total 
number of fiavours by nj = 6 and the number of colours by Nc, we also define: 

Ig = y^A - -^rii Ca = Nc. 

We denote the factorization scale by ^ip, the renormalization scale by /i/j and Euler's 
constant by "Je- The dilogarithm is defined as 

J. f . I log(l-t)dt 

Li2(2;) - ' 



t 

We split the hard matching coefficients into a representation-independent part and a part 
that depends on the irreducible representation in the colour decomposition of Appendix A. 
The nonzero hard matching coefficients for the qq — t- qq and the qq — >• qq process are the 
same provided that the appropriate representations are used. Also, as we saw in section 3, 
the contribution of the antisymmetric colour structure of the squark-squark cross section 
is suppressed for squarks that have the same flavour, so the corresponding hard matching 
coefficient is zero: 

Cqq^qq i = for cqual fiavours 
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In this equation the last two lines are proportional to the quadratic Casimir invariants of 
the representations, which are given in Eqs. (A.l) and (A. 2) for the qq — )• qq process and in 
Eqs. (A. 7) and (A. 8) for the qq — )• qq process. Furthermore we have defined the functions: 
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For the 55 — t- qq process the antisymmetric representation in Eq. (A. 4) does not 
contribute because it yields a p-wave contribution, which is suppressed near threshold. The 
hard matching coefficients for the representations of Eqs. (A. 3) and (A. 5) do contribute: 
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where in the second equation the representation / can be either of the symmetric represen- 
tations, with colour tensors given by Eq. (A. 3) or (A. 5), and the last line is proportional 
to the corresponding quadratic Casimir invariant. 

In the case of the qq — )• gg process, only the antisymmetric representation in Eq. (A. 10) 
yields a nonzero matching coefficient, the cross sections for the other representations are 
suppressed near threshold: 
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where we have defined the additional functions: 
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For the hard matching coefficients of the gg — t- gg process, the representation-dependent 
part does not scale with the quadratic Casimir invariants from Eq. (A. 11) due to contri- 
butions from box diagrams. Therefore we introduce the additional colour factors C'{Ri) 
for convenience: 
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We also introduce the function: 
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Then the hard matching coefficients are given by: 
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Note that the representation Cgg-^ggj does not contribute in SU(3). 
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The hard matching coefficients for the qg — )• qg process also have a representation- 
dependent part that is not described by a simple quadratic Casimir invariant. In fact, the 
first representation, in Eq. (A. 26), has an additional part compared to the representations 
in Eqs. (A. 27) and (A. 28), presumably due to the s-channel contribution that plays a role 
in the former. The hard matching coefficients are given by: 
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where (5/^i = 1 for the first representation and vanishes for the other representations. Also, 
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